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Quintom models, with its Equation of State being able to cross the cosmological constant 
boundary w = — 1, turns out to be attractive for phenomenological study. It can not 
only be applicable for dark energy model for current universe, but also lead to a bounce 
scenario in the early universe. 
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1. introduction 

For decades of years, the observational data has put strong evidences on existence 
of Dark Energy. The earliest evidence conies from the observation of Type la Super- 
novae by Supernova Search Team (SST) and Supernova Cosmolog y Pr oject (SCP) 
in 1998 which discovered that our universe has been accelerating \U2i. This accel- 
erating requires a kind of negative pressure matter in order to validate the Ein- 
stein Gravity. Other observations implies that our universe is nearly fiat (with 
Qtotai = 0.9996 ± 0.0199 where fltotai denotes for the total relative energy den- 
sity of our universe), while the baryon matter and cold dark matter only takes 
small part of nearly 27%, leaving large occupation for dark energy^. Due to this 
reason, it is of great importance to study the properties of dark energy. However, 
at the very beginning, people will always ask the question: what is the dark energy? 

In the literature, plenty of dark energy candidates has been proposed, see [4] 
for a review. People often classify these candidates with respect to its Equation of 
State (EoS) w = -, where p and p denotes for the pressure and energy density, 
respectively. The simplest dark energy candidate is the cosmological constant with 
energy density being near the vacuum energy p^ ss (10 _3 ey) 4 without varying 
with time^l. This candidate, proposed initially by A.Einstein, however suffers from 
the severe problem of fine-tuning and coincidence. For this sake, dynamical dark 
energy models were proposed, among which a re Q uintessence (u> > —1)^, Phan tom 
(w < —1)^, K-essence (w > — 1 or w < — 1) Quintom (w crosses — 1)^^, etc. 



1 



22,2010 19:3 WSPC/INSTRUCTION FILE proceeding 



2 Taotao Qiu 

Moreover, it is widely rea lized that accelerating can also be obtained by modifying 
the Einstein's Gravity^] However, the new released data of Supernovae, Wilkinson 
Microwave Anisotropic Probe observations (WMAP) and Sloan Digital Structure 
Survey (SDSS) as well as the forthcoming Planck etc., implies that although the 
cosmological constant with w = — 1 fits the data well, dynamical candidates still 
cannot be ruled out. Specifically, the Quintom model whose EoS can cross —1 is 
mildly favored. 

Since the observational data mildly favors Quintom, people may ask: How can 
we construct such a model theoretically? As is well known, not all of the dark en- 
ergy mode ls can make it EoS cross —1 which is constrained by the so called "No-Go 
Theorem" | 12 | 13 | 14 | 15 | 16 | 17 ] j n theorem, it is demonstrated that for theory of 
dark energy, the EoS of dark enegy model will by no means cross the cosmological 
constant boundary if it is (1) in AD classical Einstein Gravity, (2) described by sin- 
gle simple component (either perfect fluid or single scalar field with lagrangian as 
C = £((/>, d M (/>c? M </>)), and (3) coupled minimally to Gravity. Indeed, the crossing for 
such a dark energy model will lead to the divergence of physical quantities, which 
consequently result in the inconsistency of the system. For sake of this theorem, 
people proposed various kinds of Quintom models, such as double field Quintom, 
non-scalar (spinor, vector, etc.) Quintom, Quintom with higher derivative opera- 
tor, non-minimally coupled Quintom, etc. These models corresponds to violation 
of different conditions mentioned in the no-go theorem. Some of those models will 
be reviewed in detail in the next section. Due to the behavior of EoS of Quintom 
models, it can as well lead to various evolution of the universe which cannot be 
realized by non-Quintom dark energy models. 

Furthermore, if we apply the property of Quintom to the early universe, some 
interesting features will also be expected, for example, a bounce scenario will be 

II Q I 

obtained A Quintom model will give rise to a bounce scenario in AD Einstein 
Gravity, i.e. there is no need to introduce extra dimensions. This can be seen by 
investigating the conditions for a bounce to happen. In a contracting phase the 
scale factor is decreasing, i.e., a(t) < 0, while in an expanding phase the scale factor 
is increasing, a(t) > 0. Therefore we expect that at the transition point a(t) = 
while a(t) > 0. Equivalently, we expect the Hubble parameter H cross the zero 
point from H < to H > at the bouncing point, which requires the EoS of the 
universe less than —1 according to the Einstein Equation. After the bounce, in order 
for the universe to enter into the realistic one with matter dominating, radiation 
dominating, etc., it is required that the EoS of the universe being larger than — 1. 
That is, the EoS will cross —1 during the whole bounce process, and a Quintom 
behavior is needed. 

Based on the scenario of Quintom bounce, an important problem is that how can 
it give rise to the observed amount of perturbation to form the large scale structures? 
To an swer this question, we constructed the perturbation theory of Quintom bounce 
| 19 | 20| _ j n our model, Quintom field can act as inflaton after bounce and drive enough 
period of inflation. In this scenario, the inflation stage can dilute everything away 
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and set the universe to be near the state of the universe after the standard inflation, 

so the bounce process may hardly effect the following evolution of the universe and 

fits well with all the current data. Nevertheless, we do expect that bounce can 

leave some signature that can be seen for future observations. Furthurmore, using 

Quintom field we can also build another bounce model known as "Lee- Wick type" 

[9TI 

or "matter bounce", which requires no inflationary stage at all^M In this scenario, 
the perturbation is also calculated and a scale invariant power spectrum is produced 
during the contraction of the universe. 

2. the Quintom model 

Quintom model was initially proposed in |10j . where the authors combined the 
Quintessence and Phantom components together. The typical action of these kind 
of model is: 

SQmntorn = J &Xy/=g{ ^ fa Wfa ~ ^ 2 V (j) 2 - V(fa,fa)} , (1) 

where <j)\ and fa are Quintessence and Phantom components, respectively. Due to 
the combined effects of the two, it is an intuition to see that the total EoS of the 
whole system will evolve across the cosmological constant boundary. Generally, the 
effective potential can be of arbitrary form , wh ile the two components can be either 
coupled or decoupled. In the original paper the authors considered a simple form 

of decoupled case: V(fa, fa) = Vo(e m p' +e m p') where m p i denotes for the Planck 
mass while A is dimensionless constant. The evolution of the EoS was plotted in 
Fig.ffl 

After the first Quintom paper came up, lots of people investigated its property 
due to its importance. In [35], the attractor solution has been studied and in [2"5] . 
people extended to the more general case of which a coupling term has been intro- 
duced. See also [24] for more variety of double-scalar Quintom models. In [15], the 
perturbation of Quintom model has been calculated and a self-consistent perturba- 
tion theory were constructed. 

Double field is the simplest and most natural scenario of Quintom model. How- 
ever, it suffers from many problems such as big-rip and quantum instability. So 
people have to think about alternative ways to realize the EoS crossing. Another 
Quintom model is the addition of a higher derivative operator to the single scalar 
field ESI xhe most general lagrangian is as follows: 

C = C(fa V^Wfa Ufa v M v„0Vv>, ...) , (2) 

where □ = V^V^ is the d'Alembertian operator and the ellipse denotes for other 

higher dimensional operators. The higher order operators c an be deri ved from fun- 

9fi 97 9ft 

damental theories such as string theory or quantum gravity 11 1 \ and with the 
addition of high order terms to the Einstein Gravity, the theory is shown to be 
renormalizable Because of the extra degrees of freedom provided by the higher 
order term, it can simulate double-field model in some specific cases. However, it 
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has more interesting features of its own. In Higher derivative theory, the dispersion 
relation is modified, and it may provide possible solutions to the problem of quan- 
tum instability |26 | 30 ] Furthermore, because it is more complex and involved in the 
high derivative term, it may give rise to some new behaviors of evolution of the 
universe. 

We consider the phenomenological form of higher derivative action as follows 



S HD = / d A x^{A{$)V^»<t> - ^' - V(cj>)} , (3) 

J m P i 

where A(4>) and C(cp) are some functions of the field <f>. For various choice of the 
form of these functions, the evolution of the field could be very different. For sake 
of simplicity, we choose A{4>) = — | and V(<p) = as an example. By redefining the 

field variable x = ^P-Ocjj and ip = — (4> + x), we obtain a simple form of action: 

rn P i 



S eff = J ^^{--V^V^ + -V„ X V*x " c{ J +x) } > ( 4 ) 

and the last term could be viewed as an effective mass term of \ which varies in 
terms of ip and %. Here we can see, x ac ts as a normal field while ip being the "ghost" 
field with a wrong sign. We can choose the form of C((f> + x) to control the evolution 
of the two field, so that the EoS of the whole system can not only cross — 1, but 
also present novel behaviors. In our numerical analysis, we choose C(<p + x) to be 
small at the beginning and large at the end of the evolution, both of the two region 
being nearly a constant, while in the mediate region it has a significantly running. 
Therefore, at the beginning of evolution the "ghost" field ip evolves as a massless 
field with its own effective EoS nearly unity. In the mediate region, however, the 
effect of ip in potential term is involked and it behaves like a real Phantom and draw 
the total EoS below —1. In the future, it behaves like massless field again. While in 
this region the normal field x S e t s a large value of effective mass, the whole system 
will evolve as a non-relativistic matter with the EoS oscillating around zero. See 
numerical results in Fig. [2] 

For actions of field theory linear with the kinetic term, higher derivative operator 
could only exist as higher order term with an energy cut-off. However, it can also 
exist with the same order as the kinetic term in a non-linear field theory. An explicit 
example is the so-called "String Inspired Quintom" , w here such a term resides in 
the square root in the Dirac-Born-Infeld (DBI) action I^H. The effective action is 
given as: 



Sdbi = / d A x^{-V{<j } )^\ - q'V^V^ + , (5) 

where a' and /3' are coefficients of dimension 4. One can see that without the high 
derivative term, the action will be reduced to the normal DBI action describing a 
tachyon state in string theory ^1 However, one cannot remove this term by adding 
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a total integral to the action as in the linear theory. Furthermore, one can also 
include infinite numbers of higher order derivative terms as in the context of p-adic 
string theory EH. 

Due to the effect of the higher derivative operator, one can see from the numerical 
results in Fig. [3] that the EoS can also cross —1 naturally and various behaviors are 
presented according to the form of the potential. As a non-linear theory, it is also 
important to check the stability of classical perturbation. We calculated the second 
order action of this model and numerically obtained the variation of with respect 
to time during evolution. From the result one could see that for our cases c 2 s varies 
between the range of (0, 1), denoting neither instability nor unphysical propagation 
of the fluctuation. 

Of course there are other models that can make EoS crossing —1, among which 
are vector fields EI3 spinor fields ESI^ non-minimal couplin g field s E3 as well as the 
theories of modified gravity or high dimensional theories 1 ^° ! ^ which will not be 
discussed here because of the page limit. 

As a side remark, it is noticeable that due to the dynamic behavior of Quintom 
models, it can lead to many interesting fates of the universe in the future, which is 
expected to be determined by forthcoming observations. The asymptotic behavior 
of Quintom can mimic that of Quintessence, Phantom, as well as ACDM model. 
Furthermore, it can bring novel features that cannot be realized by any of them, 
for instance, the oscillating behavior around the cosmological constant bou ndary. 
Within this phenomenon, we could construct a cyclic or recurrent universe III 
Other features include that the Quintom models has a cosmic self-duality where one 
kind of Quintom model in expanding univer se is dual to the other in contracting 
universe depending on the initial conditions In this sense, Quintom model is 
of very much interest in phenomenology. Meanwhile, there are also many subtle 
issues about Quintom that remains unclear, such as its connection to the fundc- 
mental theories or particle physics, and its nature in quantum levels, etc., which are 
worthwhile of investigation in the future. 

3. Quintom bounce story: background and perturbation 

For decades of years, the theory of inflation has attracted many attentions for its suc- 
cess in solving most problems (flatness problem, horizon problem, etc) that arised in 
Standard Big-Bang Cosmology ^1. However, inflation is far fro m a complete theory 
since it suffe rs f rom other problems such as singularity problem and transplanck- 
ian problem S3. Due to this reason, people proposed several alternative solutions of 
the early universe, am ong which are Pre-Big-Bang scenario E kpyrotic scenario 



string gas scenario non-local string field theory scenario and so on. When 
restricted to 4D effective theory, a simple way to get rid of the singularity is to 
have a bou nce pr ocess at the early stage. It can be rea lized by a model of ghost 
condensate ESE3 0r the modification of Einstein Gravity ED| 51 | 52 | 53 | 54 155|. 

Quintom model, as mentioned in the introduction, can also provide a bounce 
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solution of the early universe, avoiding the singularity naturally. As an example, 
In Fig. [4] we draw a phenomenological parametrization of Quintom model. From 
the picture we can see that as the EoS evolves, the scale factor of the universe can 
transfer from damping to growing, i.e., a bounce happens. 

In order to compare the predictions of Quintom bounce scenario to the obser- 
vations, it is necessary to investigate the perturbations of the scenario. Indeed, to 
make the scenario realistic, one need the perturbations after the bounce as seeds of 
forming our galaxies and large scale structures. In singular bounce scenarios such 
as Pre Big Bang and Ekpyrosis, the fluctuation cannot evolve through the bounce 
point consistently. It will diverge at the pivot and thus invalidate the linear pertur- 
bation theory. While in our scenario, as we will see, the proper matching condition 
can be used and the fluctuations can transfer from contraction to expansion natu- 
rally. Furthermore, according to different mode, the fluctuation can possess features 
both found in singular and non-singular bouncing models. The power spectrum of 
our scenario has also been calculated, and due to different models, we can get either 
running or scale invariant power spectrum I 19 | 20 | 

In what follows, we will focus on bounce scenario caused by non-interacting 
double field Quintom. The general form of such a model is: 

£qb = ^v^v^ - iv^v^ - V(4>) - wy>) , (6) 

where V(</>) and W(ip) are potentials for normal and ghost fields. We will see that 
due to different forms of potentials, different results will come about. 
Case I. V((j>) = \m 2 (j) 2 , W{i>) = 0. 

In this case the ghost field only remains its kinetic term, thus its energy 
density evolves proportional to a" 6 and become important only near the bouncing 
point. At regions far away from the bounce, the universe is dominated by the normal 
field. We begin the contracting phase with (f> oscillating around the minimum of the 
potential. Because of the contraction of the universe, the amplitude of the oscillation 
grows as a 3 / 2 and the field behaves as non-realistic matter. This period is called 
"Heating phase". After the last oscillation, the field climb up along the potential 
and caused a period of "deflation" of which the EoS of the universe approximately 
near —1. Meanwhile, the energy density of the field if) is growing all the time. When 
it catches up with that of the field 0, the total energy density vanishes. As we learn 
from Einstein Equation that H = with a positive time derivative H . Thus the 
bounce happens. After the bounce the ip field damps quickly while the (f> field rolls 
down slowly along its potential, very much like the chaotic inflation model. At last, 
4> oscillates around the minimum of the potential again, with a damping amplitude. 

The left hand side of Fig.[5]is the evolution of EoS during the bounce process and 
the right hand side is the sketch plot of the space-time in bounce scenario while the 
horizonal and vertical axis denotes for physical distance and time respectively. The 
black solid line represents the Hubble horizon and green and blue line are different 
fluctuation modes. 
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From perturbed Einstein Equation, we get the equation for the Newtonian grav- 
itational potential $ as 

$" + 2{H - Q-W + 2{H' - U^-)<5> - V 2 $ = 8ttG(2H + 4r) W , (7) 

(jy (fr cfy 

where % is the comoving Hubble constant, and prime denotes derivative with respect 
to comoving time rj. For the initial condition, we set it as Bunch-Davies vacuum 
as usual in the far past, when u = ^ ~ / „ and $ ~ ri~ 3 . 1 „ . Here we use the 

1 ' <t>' V2k* 1 V2k? 

fact that at the very beginning when w«0 we have — ~ r]~ 3 . At regions far from 
bouncing point, the right hand side of Eq. (J7J can be neglected and the equation 
becomes homogeneous. It has two branches of solution at super-Hubble region, one 
is growing while the other is constant. At the regions near the bounce, we get both 
of the branches oscillating with the amplitude depending on the energy scale of 
the bouncing. After the bounce, we again get two branches where one is constant 
while the other is damping. The solution s of ea ch stage is matched consistently via 
Deruelle-Mukhanov matching conditions 1 5 1 6 1_ 

Unfortunately, this kind of model cannot give rise to a scale invariant power 
spectrum as obtained by observation. The reason is easy to see: In the contracting 
phase, the perturbation will get out of the horizon and the amplitude of the per- 
turbation will be raised by the growing mode. Such a behavior corresponds to the 
modification of the initial condition of perturbations at the following inflationary 
stage. The power spectrum will thus have a blue tilt. For the above reason, some 
alternative models need to be considered. 

Case II: V{<t>) = M - \) + ^, W(tl>) = 0. 

In this case^^ni, the potential of normal field <p is of Coleman- Weinberg type, with 
two vacua set on both sides of the middle see Fig. [6] We set our initial state with 
4> residing in one of the vacua while oscillating. When the amplitude grows large 
enough the field climbs to the plateau of the potential, while the growing energy 
density of the second field ip catches up, and the bounce happens, followed by a 
period of inflation. As one could see from the numerical analysis of background 
evolution that, there is no "deflation" period in the contracting phase and the 
symmetry with respect to the bouncing point has been broken, see Fig. [7j Thanks 
to this symmetry violation, we could see from the sketch plot that for large k modes, 
the fluctuation will stay inside the horizon all the contracting time, leading to a scale 
invariant power spectrum. Meanwhile, for small k modes which exits horizon during 
contraction, the spectrum will have a blue tilt. 

Case III: V{4>) = ±to 2 </> 2 , W(ip) = -±MV- 

This is another interesting case where both of the field get a mass term, while 
that of the ghost field has a wrong sign 1 ^. This model can be obtained from the 
scalar sector of Lee- Wick Standard Model with M ^> m ^21. In this case, both of the 
fields start off oscillating around their extremes of the potential, and the initial stage 
is dominated by the normal field. The oscillation amplitude of both fields scale as 
a(t)~ 3/2 and the universe presents non-relativistic-matter-like. However, ip oscillates 
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much rapidly due to the heavy mass, and eventually the energy density of tp catches 
up with <j) and bounce happens. After the bounce, both fields continues oscillating 
with a damping amplitude. During the whole process except the bouncing point, 
the EoS of the universe has an average value of w = and there is neither deflation 
nor inflation era at all. See Fig. [H 

In this model, a scale invariant power spectrum can be obtained during contract- 
ing phase for fluctuations on scales larger than Hubble radius. This could be seen 
easily if we work in terms of the como ving cu rvature perturbation or equivalently, 
the Mukhanov- Sasaki variable v = | 63 | 64 | w ]^ ere z ^ a for time- independent EoS. 
The well-known equation of motion for v is: 

v" + (k 2 -—)v = 0. (8) 

z 

For super-Hubble fluctuation modes, one can neglect the k 2 term and make use 
of the parametrization a(rj) ~ rf for w — to get v ~ rf 1 , and hence obtain the 
fc-dependence of power spectrum: 

p c (k, n) ~ k 3 \v(v)\ 2 a(v)- 2 ~ fc 3 K^(fc))| 2 (^) 2 ~ fc 3 - 1 - 2 ~ fc° . (9) 

V 

In deriving the above formula we've used the fact that the comoving time when 
fluctuations cross the Hubble horizon Tjff(fe) ~ fc -1 . 

4. summary 

This talk is mainly focused on the properties of Quintom dark energy model and 
its application to the early universe. Quintom model, which needs multi degrees 
of freedom, has an EoS crossing — 1 and can bring various features of the universe 
in the future. When applying to the early universe, it can give rise to a bouncing 
scenario. The perturbation theory of Quintom bounce is self consistent and a scale 
invariant power spectrum can be obtained. 
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Fig. 1. The evolution of the effective EoS of double- field Quintom model proposed in [10]. The 
parameters are chosen as: Vq = 8.38 X 10~ 12e m^, A = 20. The initial conditions were chosen as: 
4>u = -1.7m pl , <f>2i = -0.2292m pi . 





Fig. 2. The evolution of EoS of the model proposed in [31]. For the left hand side case the initial 
values are i/>i = -0.26m pJ , ipi = 3.52 X lO" 62 ™^, X i = 0.25m p! , \i = -3.62 X lO -62 ?™ 2 ,, while 
for the right hand side case they are ipi = — 0.26m p j, ipi = —2.84 X 10~ 62 m 2 ; , \i = 0.25m p ;, 
Xi = 2.74 X 10- 62 m 2 . 
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Fig. 3. The first row: The evolution of the EoS of the "String Inspired Quintom model" , of which 
the first picture for potential V(<j>) = Vbe — ^ while the last two for potential V(<j>) = e \</,^°-\ < t > 
and different parameter choices. The second row: The variation of the corresponding c? s of each 
case obove. 




Fig. 4. The evolution of the EoS w, hubble parameter H and the scale factor a as a function of 
the cosmic time t. The EoS is parameterized as w = —r — ^ where r = 0.6 and s = 1. 




Fig. 5. The left hand side: The evolution of the EoS in Quintom bounce model of Case I. The 
initial values are chosen as <f>i = —5.6 X 10~ 3 m pi , <j>i = 2.56 X 10~ 10 nip l , ipi = 4.62 X lO -85 ™.^ 
with m = 1.414 X 10 — 7 m p ;. The right hand side: A sketch of the evolution of perturbations with 
different comoving wave numbers k in this case. 
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Fig. 6. The sketch plot of the potential as function of <f> in Quintom bounce model of Case II. 




Fig. 7. The left hand side: The evolution of the EoS in Quintom bounce model of Case II. The 
initial values are chosen as tj>i = — 0.82m p ;, ifc = 3.0 X lO -10 ™,^, ipi = 5.0 X W~ 13 m 2 ll with 
A = 8.0 X 10 — 14 and v = 0.82m p ;. The right hand side: A sketch of the evolution of perturbations 
with different comoving wave numbers k in this case. 




Fig. 8. The left hand side: The evolution of the EoS in Quintom bounce model of Case III. The 
initial values are chosen as fa = 1.74 X 10~" 3 m pI , (f>i = 1.44 X lO" 8 ?™^, ipi = 8.98 X 10 _6 m p ;, 

^ = -14.08 X lQ- 12 m 2 pl with m = 5.0 X 10~ 6 m pi , M = 1.0 X 10~ 5 m pi . The right hand side: A 
sketch of the evolution of perturbations with different comoving wave numbers k in this case. 



